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We show that a Fermi surface in underdoped YBa2Cu306+ x yielding the distribution of quan- 
tum oscillation frequencies observed over a broad range of magnetic field can be reconciled with 
the wavevectors of charge modulations found in nuclear magnetic resonance and x-ray scattering 
experiments within a model of biaxial charge ordering occurring in a bilayer CuC>2 planar system. 
Bilayer coupling introduces the possibility of different period modulations and quantum oscillation 
frequencies corresponding to each of the bonding and antibonding bands, which can be reconciled 
with recent experimental observations. 
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Our level of understanding of metals typically depends 
on the extent to which Fermi surface cross-sections ob- 
served in quantum oscillation experiments can be ex- 
plained by comparison with band structure calculations, 
and the degree to which Fermi liquid behavior prevails [Tj- 
4 . Additional factors, such as a periodic modulation of 
the electronic states — caused, for example, by a spin- or 
charge- density wave — can introduce additional gaps, in 
which case the Fermi surface needs to be understood in 
terms of the original band structure translated by order- 
ing wavevectors @HZ]- 

In the case of high temperature superconductors, a de- 
termination of the electronic structure is essential for un- 
derstanding the origin of pairing. Yet, the small Fermi 
surface pocket observed in the underdoped cuprates 
YBa 2 Cu30 6 +a; [5J and YBa 2 Cu40 8 [H] cannot readily 
be explained in terms of the unmodified band struc- 
ture [5Ul7j. The momentum-space cross-section of the 
experimentally observed pocket is roughly thirty times 
smaller than the size of the large hole Fermi surface sec- 
tions expected for the Cu02 planes from band structure 
calculations 18-20J, with the cyclotron motion having 
apparently the opposite sense of rotation to that ex- 
pected [10l [12]. A small pocket has been proposed to 
originate from band structure involving the mixing of the 
chain states with BaO orbitals in YBa2Cu 3 06+ s [Ti?ll20) . 
but the relevant band is reported to lie significantly 
(sa 0.6 eV) below the Fermi energy in photoemission ex- 
periments [TT] . 

The emergence of long range charge order reported 
in underdoped YBa2Cu306+ x by nuclear magnetic res- 
onance (NMR) [21] and x-ray scattering |22"24] may 
contribute to the observed difference between the small 
Fermi surface pocket measured and the large Fermi sur- 
face expected from band structure. A direct link, how- 
ever, has yet to be established between the structure 
of the measured order and reconstruction of the large 
Fermi surface to yield the observed small pocket. Fur- 
thermore, some differences are reported in the ordering 



wavevectors seen in different experiments [2D - f2"4] . and it 
remains to be settled whether the observed charge or- 
der in YBa 2 Cu 3 6+a; is stripe- like [2"5H2"7] as opposed to 
biaxial [T7] [55] in nature. 

A satisfactory electronic structure model of under- 
doped YBa2Cu30g+a; would be expected to explain the 
ordering wavevectors found in both NMR and x-ray 
experiments, while accounting for the measured distri- 
bution of bilayer-split frequencies in quantum oscilla- 
tion experiments [111 as well as other experimen- 
tal observations relating to the electronic structure [TD1 - 
H3 ED 113 [SI [S3]- A single layer scheme of Fermi sur- 
face reconstruction by biaxial charge order has previously 
been shown to yield nodal electron pockets [15 ] fT7 ] l28 | [35] 
consistent with experimental observations of the small 
size of the electronic heat capacity at high magnetic 
fields [301 131] > the negative electrical transport coeffi- 
cients [TUHH] and the predominantly nodal Fermi surface 
location inferred from photoemission, Fermi velocity and 
chemical potential oscillation measurements [HI [17l [33] . 
In this paper, we show that the introduction of bilayer 
coupling [33] to the single layer scheme [TH [T71 [251 [35] 
leads to the possibility of a modulation period and Fermi 
surface pocket of different size corresponding to each of 
the bonding and antibonding bands, consistent both with 
experimental measurements of the charge ordering |21| - 
|2"4"] , and quantum oscillations [111 [2S] ■ 

NMR experiments provide evidence for long range 
charge ordering in a magnetic field B gl 15 T [3T] in 
underdoped YBa2Cu306+a;, while recent x-ray scatter- 
ing experiments reveal ordering vectors that are two- 
dimensional [22 24 . The x-ray results are significant for 
two reasons. First, they reveal the scattering intensity 
and correlation lengths to be equivalent along the a and 
6 lattice directions of the orthorhombic unit cell [T5] , with 
the orthorhombicity and CuO chain oxygen ordering hav- 
ing no apparent effect in causing stripe domains to form 
preferentially along either a or b. Second, they reveal the 
charge ordering to be incommensurate with wavevectors 
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Qa = (j^, 0, Qz) and Q b = (0, jfr, Q z ) [of in-plane pe- 
riod A a = Afc w 3.2 lattice constants] consistent with a 
wavevector nesting the bilayer-split bonding band [22H24] 
(where Q z refers to the c-axis component of the ordering 
vector). The periods reported in x-ray scattering exper- 
iments, however, differ from that (A a — 4.0) reported in 
NMR studies [5T] of similar YBa2Cu306+:r samples and 
that (A « 4.2) [SS] typically found by x-ray diffraction in 
the same region of the phase diagram in single layer spin- 
stripe ordered cuprates [37] . Furthermore, no obvious 
relationship is found between the charge ordering wave 
vectors found in x-ray experiments on YBa2Cu30g+x and 
low energy spin excitation wave vectors |38j . 

Given the unambiguous planar origin of the two- 
dimensional charge ordering found in x-ray scattering 
experiments in YBa2Cu306+a; [24] . we focus here on a 
model dispersion solely for the CuC>2 bilayers. While 
photoemission experiments have yet to identify signa- 
tures of band folding consistent with charge ordering, the 
wave vectors found using x-rays in YBa2Cu3 06+a; (22] 
153] . as well as those found earlier in doping-dependent 
scanning tunneling microscopy measurements in other 
cuprates [3U 20], suggest Fermi surface nesting at the 
antinodes [35] [39 . The opening of an antinodal gap 
(shaded light grey in Fig. [TJ that is expected to accom- 
pany ordering at such wave vectors is consistent with a re- 
constructed Fermi surface that involves relative transla- 
tions of the remaining arc-like segments of Fermi surface 
(represented by thick black lines in Fig. [T] [T71 12"51 135] ). 
These translated segments of Fermi surface would oc- 
cupy the same regions in momentum space where 'Fermi 
arcs' are observed in photoemission experiments in un- 
derdoped cuprates [HJ |4"0"H4"2"] . 

The single layer unreconstructed Fermi surface pro- 
duced using e(k) = £o + 2ti [cos ak x + cosbk y ] + 
2tn[cos(ak x + bk y ) + cos(ak x — bk y )\ + 2t 2 o [cos 2ak x + 
cos 2bk y ] has a single nesting wavevector spanning the 
antinodal regions [i.e. |Q Q | = |Qb| in Fig.^a)]. On intro- 
ducing bilayer splitting, two nesting vectors now span the 
bonding and antibonding bands such that |Q B | = _|Q B \ 



and |Qa B | = |Q^ B | are inequivalent [shown in Fig.[2jb)]. 
For the single-layer dispersion [43], we choose t\\/tw = 
— 0.32 and ^o/^io = 0.16 after Andersen et al. [TH] and 
Millis & Norman [25]. The bonding and antibonding 
bands are obtained by diagonalizing 
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after Garcia- Aldea and Chakravarty [34], where 
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and t c represent intrabilayer and interbilayer couplings 
respectively. On evaluating Eq. ([!]) one obtains 

e(k) B ' AB = e(k) t V^( k ) 2 +t*+2t ± t c coack z (2) 




FIG. 1: Single layer schematic showing how gapping of the 
electronic density of states at the antinodes (depicted in grey) 
leaves remaining arc-like segments of Fermi surface [1711281135] 
(depicted in black), whose reconnection upon translation by 
Q a and Qb leads to the possibility of a nodal pocket for which 
the sense of cyclotron motion in a magnetic field is opposite 
to that of the original hole sheets. Thick black lines represent 
e(k), thick magenta lines represent e(k + Q a ) and e(k + Qb), 
while thick cyan lines represent e(k + Q a + Qb). No direct 
coupling occurs between e(k) and e(k+Q a + Q(,) in the model, 
however, e(k) and e(k + Q a + Qt) are both coupled to e(k + 
Q a ) and e(k + Qb), thus leading to a closed pocket. 
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FIG. 2: (a) Single layer Fermi surface corresponding to e(k) 
as defined in the text, (b) Bilayer Fermi surface after Eq. (J2| 
for k z = 0, where the red and blue lines correspond to bond- 
ing and antibonding bands respectively. The cross-sectional 
area of the Fermi surface represents an effective hole doping 
of p = 8 % relative to the half filled band, which we estimate 
from the total width of the chain Fermi surface seen in photoe- 
mission experiments on ortho-II ordered YBa2Cu306+a; [llj . 
Q B and Q AB refer to nesting wavevectors of the bonding 
and antibonding bands respectively, while Q a and Q;, per- 
tain to the different lattice directions in the orthorhombic 
unit cell |18j . £ gap refers to the nodal energy gap. 



for the bonding and antibonding bands (shown for 
cos ck z = in Fig. |2|d) , where c refers to the c-axis lat- 
tice parameter. Setting tj_oAio ~ 0.5 and t c /t±o = 0.25 
produces a bilayer splitting at the antinodes [k ant i n odc = 
(±^,0, k z ) and (0, k z )] resembling that in bandstruc- 
ture calculations [HI [20] in YBa2Cu306+a;, but with a 
reduced nodal gap (e g ~ 2t c ) that falls within the exper- 
imental upper bound from photoemission [42j and mag- 
netic quantum oscillation measurements 16J. 
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FIG. 3: Schematic of the interlayer Fermi surface corruga- 
tion at bk y = ±7r according to Eqn Q. Owing to a Q z = 7r 
component, translation (dotted lines) of the bonding and an- 
tibonding Fermi surfaces by Q B and Q AB , respectively, leads 
to almost perfect interlayer Fermi surface nesting at the antin- 
odes. 



The corrugation of the unreconstructed bilayer Fermi 
surface with respect to the interlayer momentum k z 
(shown in Fig. [3]) caused by t±o and t c implies that 
Fermi surface nesting is optimized by acquiring a Q z = tt 
component to the wave vector. Such a component pro- 
duces a reconstructed Fermi surface in which the fun- 
damental component of the corrugation is suppressed, 
leaving a small residual second harmonic corrugation 
of order t^./t±Q [originating from the square root in 
Eqn. Q]. Hard x-ray scattering experiments find di- 
rect evidence for an interlayer component to the order- 
ing vector |12j . Furthermore, quantum oscillation exper- 
iments performed over a broad range of magnetic field 
and angles are found to be consistent with a bilayer-split 
Fermi surface that exhibits a comparatively weak corru- 
gation |15[ 116] . The strong suppression of the interlayer 
optical conductivity within the underdoped regime [U] 



could be a further signature of Fermi surface corruga- 
tion suppressed by nesting. The small size of the resid- 
ual corrugation term (which we estimate later in the 
manuscript) compared to the other energy scales enables 
us to consider a simplified two-dimensional approxima- 
tion to the Fermi surface, which we obtain by neglecting 
the fcz-dependence in Eqn 

To obtain the full set of reconstructed bands, the rel- 
evant couplings (V^ AB ) between bands with relative 
translations by ±Q B , ±Q B , ±Q AB and ±Q^ B need to 
be considered [including both the bonding and antibond- 
ing bands given by Eqn ^2j], requiring the construc- 
tion of a Hamiltonian in the form of a matrix (see e.g. 
Refs. [251131]). Orthogonality between the Bloch wave- 
functions of the bonding and antibonding bands ensures 
that charge-density modulations form independently for 
those two bands, with the coupling between them being a 
second order effect [45]. The strongest coupling V ,6(k') 
is therefore expected to occur between states within each 
of the bands translated by wavevectors intrinsic to that 
band — i.e. between e(k) B and e(k±Q B h ) B and between 
e(k) AB ande(k±Q AB ) AB 

While the complete set of bands for strictly incommen- 
surate order (in which A B and A AB are irrational) requires 
an infinite number of translations and a matrix of infi- 
nite rank, which cannot be diagonalized, it was shown 
in Ref. [3 5) that one can approximate incommensurate 
behavior by considering rational values of A (in which A 
is the ratio of two integers) — yielding large but manage- 
able matrices. Provided <C |V 0) b(k')| < |tio|, the Fermi 
surface consists of multiple repetitions of the same elec- 
tron pocket throughout the Brillouin zone [35] ■ Here we 
go a step further in simplification by considering only the 
minimal number of terms necessary to produce a single 
instance of the electron pocket for each band (e.g. those 
illustrated in Fig. [I]). Hence, 
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which must be diagonalized for both bonding and anti- 
bonding bands. With these simplifications, we can esti- 
mate the size of the pockets, their effective masses and 
contributions to the electronic density of states for arbi- 
trary values of A B or A AB — bearing in mind that the 
complete Fermi surface will have these same pockets re- 
peated by ±Q B , ±Q B , ±Q AB and ±Q AB and combina- 



tions thereof throughout the Brillouin zone. 

We procede to obtain the Fermi surface for a given 
strength (Vq) of the potential by determining the min- 
imum in the electronic density of states (DOS) for the 
bonding and antibonding bands as a function of Q B b and 
Qaf ■ We utilize the fact that the DOS is further mini- 
mized by considering momentum-dependent couplings of 
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FIG. 4: Bonding (red) and antibonding (blue) nodal electron 
Fermi surface pockets obtained using Eq. |3| for Vo/t — 0.7, 
depicted in the extended Brillouin zone representation. Black 
lines indicate the unreconstructed Fermi surface shown in 
Fig. [2jb). Grey lines represent the remaining unnested por- 
tions of Fermi surface resulting from the inclusion of only the 
lowest order terms in Eq. (|3|. The grey lines become recon- 
structed into multiple repetitions of the pockets throughout 
the extended Brillouin zone on including more terms in the 
Hamiltonian [5S]. A schematic representation of the com- 
plete reconstructed Fermi surface is produced by translating 
the closed pockets by ±Q B , ±Q B , ±Qa B and ±Q^ B and all 
multiples and combinations thereof (not shown for clarity). 
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FIG. 5: (a) The solid lines (left-hand axis) show the optimal 
values of A B and A AB at which the electronic density-of-states 
(DOS) is minimized, which are found to depend on the cou- 
pling strength Vo. Both the bonding (B) and antibonding 
(AB) bands are shown. At Vo/iio = 0.7, the values of A agree 
with those obtained from x-ray scattering [22] and NMR [21] 
experiments (filled circles). The dotted line (right-hand axis) 
shows the corresponding pocket frequencies calculated for the 
bonding and antibonding bands [related to the pocket area 
via Onsager's relation F — (h/2ne)Ak], which also depend on 
Vo- At Vb/tio = 0.7, the frequencies coincide with the bilayer- 
split extremal orbits Fq ± 2AF (filled squares) inferred from 
quantum oscillations in which the frequencies Fo and Jb± AF 
are observed JTB]. (b), Electronic DOS for Vo/t w = 0.7 (plot- 
ted in units of the band mass mb per band in free electron 
masses m c ) at the Fermi energy as a function of 1/A, reveal- 
ing a clear gap for each band with a minimum that is used to 
infer the optimal value of A. Red and blue lines correspond 
to the DOS of the reconstructed bonding and antibonding 
bands obtained from Eq. J3J|. Dots indicate the minimum 
DOS per band within the reconstructed Brillouin zone when 
all multiples of Q B ',f B are included in Eq. ^ [46] . 



the form [35] 

V a (k) B < AB = V 1=^(1 - r B < AB cos bky) 
Vb{k ) B ' AB = V 0J ^ w (l-r B > AB cosak x ). (4) 

In real-space, r corresponds to a modulation of bond 
strength, for which there is evidence in recent experi- 
ments [23]. The value of r is tuned so as to suppress the 
energetically unfavourable coupling between translated 
states whose group velocities point in the same direction. 
By adjusting Vq, tj_ , r B and r AB (where B and AB refer 
to the bonding and antibonding bands) , we are then able 
to tune the locations in A B and A AB of the minima in the 
DOS. 

We find that the wavevectors Q B fc and Q AB at which 
the electronic DOS is minimized when Vo/iio = 0.7, 
i-LoAio = 0.54, r B ps 1.65 and r AB pa 1.35 correspond 
closely to the wavevectors reported in x-ray scattering 
and NMR experiments. These same wavevectors also 
produce nodal pockets (see Fig [I| with areas consis- 
tent with those found in recent quantum oscillation ex- 
periments uncovering a bilayer-split Fermi surface [see 
Fig. p)Fa)]. Figure [5jb) shows the electronic DOS (solid 
lines) for Vb/iio = 0.7 plotted as a function of 1/A 



for both the bonding and antibonding bands in units 
of the effective mass. The deep minima (where the 
free energy is minimized) are located at X B b = 3.30 
and A AB = 4.05, which are very similar to the values 
X a ,b = 3.28 and A a = 4.00 reported in x-ray scatter- 
ing [22H23] and NMR [5T] experiments, respectively, per- 
formed on underdoped YBa2Cu306+ic [46] . 

Using Onsager's relation, F = (H/2ne)Ak (where Ak 
is the Fermi surface orbit cross-section in momentum 
space), we obtain F B ps 337 T and F AB ps 709 T 
quantum oscillation frequencies for the model Fermi sur- 
face pockets in Fig. [4] when Vb/tio = 0.7 (see Fig. |5j. 
The size of the closed orbits is not found to be signifi- 
cantly impacted by the inclusion of higher order terms 
in the Hamiltonian in the limit |V a ,&(k')| < t 10 , justi- 
fying our neglect of such terms for estimating the orbit 
areas and their effective masses. The observed frequen- 
cies would be further modified due to magnetic break- 
down tunnelling, which is expected due to the suppressed 
gap e gap pa 10 meV between bonding and antibonding 
bands at the nodes in underdoped YBa 2 Cu30 6+a; , as 
seen in photoemission experiments [42 . The low ef- 
fective mass of mo ~ 0.34 m G for the prominent fre- 
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quency To from band structure calculations [151 |2"5] due 
to the hopping parameter tin = 380 meV implies a 
renormalization of the bands (and subsequently Vo) by 
a factor of w 4.7 in the vicinity of the Fermi energy 
for consistency with the experimentally measured value 
77i* w 1.6 m c . Such a value of the renormalization 
yields for the residual corrugation (neglected in the two- 
dimensional Fermi surface approximation) (t;?/t_i_n)/4.7 
ss 2.5 meV, which is comparable in magnitude to the 
small interlayer coherence temperature (~ 27 K) reported 
in interlayer transport experiments |47j . Finally, we ob- 
tain V = 0.7 x 380/4.7 « 57 meV for the strength of the 
charge order parameter, which, on taking To f=a 135 K [23] 
for the ordering temperature, yields 2Vq/ ksTg ss 10. 
While this ratio is large, it lies within the range of val- 
ues, 6 < 2Vo/ TbTo < 17, found exprimentally in typical 
charge- density wave systems [4"8H52"] . 

The magnetic breakdown resulting from the small gap 
£gap between bonding and antibonding bands at the 
nodes (see Fig. [6]) is expected to occur at four points 
around each pocket (in the repeated Brillouin zone), giv- 
ing rise to a sequence of five frequencies Tn+mAT (where 
?7i = —2, —1, 0, 1 and 2) in which the original bilayer- 
split frequencies correspond to T B = T — 2AT and 
T AB = T + 2 AT [IB] (see Fig. A small magnetic 
breakdown gap (e ga p) would further cause the central 
frequency Fq (depicted in Fig. [6]) to dominate [IB] with 
weaker adjacent features at To — AT and To + AT [16], 
while the spectral features at T ± 2AT are expected to 
be further suppressed in amplitude [IB]. The calculated 
values To ss 519 T and AT ss 95 T compare favorably 
with those T « 532 T and AT w 90 T found experi- 
mentally [16 . 

We turn to the question of why the wavevectors for 
the bonding and antibonding bands might be detected in 
different experiments [51] [55] and consider whether de- 
tails of the Fermi surface reconstruction could provide 
potential clues. X-ray scattering experiments are ex- 
pected to be sensitive to the momentum-dependent sus- 
ceptibility [55H5I] , which is largest for the sections of the 
Fermi surface that are better nested and most completely 
gapped once Fermi surface reconstruction takes place. 
We find the electronic DOS to be lower for the bonding 
band (reflecting its flatter topology) , which could poten- 
tially account for the observation of A B in x-ray scattering 
experiments [22-24.. Further experiments are required 
to determine whether a secondary feature can be seen 
with x-rays corresponding to the antibonding band at 
A AB w 4; possibly induced by the application of a mag- 
netic field at low temperatures [53] . NMR, by contrast, 
is expected to yield a splitting when the charge modu- 
lations become commensurate with (or locks-in to) the 
underlying crystalline lattice [21] ; owing to its proximity 
to an integer value, such a possibility is more likely for 
A AB . 

We also note that whereas the formation of a Cu02 




FIG. 6: (a) Schematic showing how magnetic breakdown 
tunneling (black dotted line) leads to the observation of 
the To frequency, which dominates the quantum oscillation 
Fourier spectrum, considering the hypothetical case where 
the bonding and antibonding charge ordering vectors are 
the same (i.e. Q B = Q AB ) after Ref. [16]. On the left- 
hand-side, magnet breakdown tunneling of probability p 2 = 
exp(—e 2 B/h 2 uj 2 F) [2] (where cj c = eB/m* is the cyclotron 
frequency) across the nodal gap e gap separating bonding and 
antibonding Fermi surfaces (indicated in red and blue respec- 
tively) is seen to occur at four equivalent points by symmetry 
around each of the reconstructed pockets. The resulting mag- 
netic breakdown orbit is indicated on the right-hand-side, (b) 
On considering the case relevant to the present paper in which 
Q B 7^ Q AB , several pockets with relative momentum space 
translations of ±(Q B - Q AB ) are incorporated into the mag- 
netic breakdown path (left-hand-side schematic). For clarity, 
only those pockets directly involved in magnetic breakdown 
tunneling are shown (the complete Fermi surface would con- 
sist of the repeated translation of these pockets by all possi- 
ble combinations of Q B , Q B , Q AB and Q AB ). The magnetic 
breakdown orbit arising from the four equivalent symmetry- 
related points of magnetic breakdown tunneling (neglecting 
smaller gaps) in this case has the same frequency To (right- 
hand-side) as that in (a). 



plane-centered charge modulation involving the bond- 
ing band [23] is mostly decoupled from oxygen ordering 
within the chains, as indicated by recent x-ray scattering 
experiments |24| . the opposite may be true for a charge 
modulation involving the antibonding band owing to a 
greater overlap of its wave function with the chains (the 
chains themselves being prone to different forms of or- 
der [53] [55J). Local interactions with the oxygen order 
could cause the latter to become more susceptible to a 
lock-in transition or to differences in the strength of the 
modulation along the a and b lattice directions, poten- 
tially revealed in NMR experiments [21) . 

In summary, we have shown that the different wavevec- 
tors reported in x-ray scattering and NMR experiments 
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and the frequencies observed in quantum oscillation ex- 
periments [16] can be consistently explained by biaxial 
charge ordering involving different nesting vectors for the 
bonding and antibonding bands in the bilayer cuprate 
YBa 2 Cu 3 6+a; PU HH [M] ■ We note that the negative 
electrical transport coefficients, size of the high magnetic 
field electronic heat capacity and chemical potential os- 
cillations are similar to those expected in previously con- 
sidered biaxial charge ordering models in which the ef- 
fects of bilayer coupling were neglected [TS1 HH1 ES] • The 
high purity of YBa2Cu3C>6+:r samples together with their 
suitability for both quantum oscillation and x-ray scat- 
tering experiments make these materials a model system 
for understanding the interplay between charge ordering 
and superconductivity in the cuprates. 
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